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A systematic fieldtheoretic description of critical adsorption on curved 
objects such as spherical or rodlike colloidal particles immersed in a fluid 
near criticality is presented. The temperature dependence of the corre- 
sponding order parameter profiles and of the excess adsorption are calcu- 
lated explicitly. Critical adsorption on elongated rods is substantially more 
pronounced than on spherical particles. It turns out that, within the con- 
text of critical phenomena in confined geometries, critical adsorption on 
a microscopically thin 'needle' represents a distinct universality class of its 
own. Under favorable conditions the results are relevant for the flocculation 
of colloidal particles. 
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I. INTRODUCTION 



In colloidal suspensions the interaction between the mesoscopic dissolved particles 
and the solvent is of basic importance [MB- F° r example, the solvent generates effective 
interactions between the colloidal particles which can even lead to flocculation. The 
richness of the physical properties of these systems is mainly based on the possibility to 
tune these effective interactions over wide ranges of strength and form of the interaction 
potential. Traditionally this tuning is accomplished by changing the chemical composition 
of the solvent, e.g., by adding salt, polymers, or other components Q. Compared with 
such modifications, changes of the temperature or pressure typically result only in minor 
changes of the effective interactions. This, however, is only true as long as the solvent is 
not thermodynamically close to a phase transition of its own. For example, if the solvent 
consists of a binary liquid mixture close to a first-order demixing transition into a A-rich 
and a B-rich liquid phase, even slight changes of the temperature or of the partial pressures 
of the two species A and B can lead to massive changes of the effective interactions between 
dissolved colloid particles induced by the occurrence of wetting transitions. They lead to 
wetting films of the preferred phase coating the colloidal particles ||. These wetting films 
can snap into bridges if the particles come close to each other leading to flocculation 
HH. For charged colloidal particles such as silica spheres immersed in the binary liquid 
mixture of water and 2,6-lutidine [Q flocculation can also be influenced by screening effects 
generated by the adsorbed layers 0. 

Similarly drastic effects can occur if the solvent is brought close to a critical point. 
The inevitable preference of the surfaces of the colloidal particles for one of the two 
solvent species of a binary liquid mixture near its critical demixing point or for the liquid 
phase of a one- or two-component solvent fluid near its liquid- vapor critical point results 
into the presence of effective surface fields leading to pronounced adsorption profiles of the 
preferred component. This so-called 'critical adsorption' becomes particularly long-ranged 
due to the correlation effects induced by the critical fluctuations of the order parameter 
of the solvent. In the case of a planar wall critical adsorption has been studied in much 
detail [fT[-|T6|j. Asymptotically close to the critical point T c it is characterized by a surface 
field which is infinitely large so that the order parameter profile actually diverges upon 
approaching the wall up to atomic distances a. As compared to a planar surface critical 
adsorption on a spherical particle is expected to exhibit important differences in behavior 
because the confining surface has a positive curvature and because a sphere represents only 
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a quasi-zero-dimensional defect floating in the critical fluid. The interference of critical 
adsorption on neighboring spheres gives rise to the so-called critical Casimir forces [17] 
which have been argued to contribute to the occurrence of flocculation near T c [|l^-|20H- 
A quantitative understanding of these phenomena requires the knowledge of the critical 
adsorption profiles near the colloidal particles and the resulting effective free energy of 
interaction in the whole vicinity of the critical point, i.e., as function of both the reduced 
temperature t = (T — T c )/T c and the field h conjugate to the order parameter. This 
ambitious goal has not yet been accomplished. Instead, the introduction of a surface 
curvature has limited the knowledge of the corresponding critical adsorption so far to 
the case of spheres for the particular thermodynamic state (t, h) = (0, 0) of the solvent 



21.22]. Only recently at least the temperature dependence of the critical Casimir force 



between a sphere and a planar container wall has been addressed ] 23f| . Thus the present 
study of the temperature dependence of critical adsorption on a single sphere contributes 
one step towards reaching the aforementioned general goal. 

Apart from spherical particles, also rodlike particles play an important role |24| . Rod- 
like objects are provided, e.g., by fibers or colloidal rods ||24|] , semiflexible polymers with 
a large persistence length such as actin [25], microtubuli f2"5 |, and carbon nanotubes ||26|| . 
Moreover the knowledge of the general curvature dependence of critical adsorption is also 
relevant, e.g., for curved membranes [0,0 dissolved in a fluid near criticality or for the 
liquid-vapor interface between a binary liquid mixture near its critical demixing point 
and its noncritical vapor, which exhibits rippled configurations due to the occurrence of 



capillary waves ]27J . 

Near criticality the relevant length scales of the solvent structures are dominated by 
the diverging bulk correlation length ^± = ^\t\~ u , where v is the standard universal 
bulk critical exponent; £g~ and £j are nonuniversal amplitudes in the one- (+) and two- 
phase region (— ), respectively, with values typically in the order of a few A. In practice 
the correlation length can span the range between 5 A - 1 /im depending on t. In the 
present context this length scale is played off against the length scale R of the radius 
of the dissolved particles. We note that the available systems can realize both the limit 
R/£ > 1 as well as the opposite limit R/£ <C 1. In the case of Ludox silica particles 
R « 12 nm |2£| so that the limit R/£ <C 1 can be easily achieved even with the upper 
limits for £ set by finite experimental resolutions. The ratio of the length I and the 
radius R of rodlike particles can be quite large, in conjunction with a small radius such 
as R ~ 7nm in the case of colloidal boehmite rods pi. In this work we consider long 
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rods, i.e., R, £ -C I, and neglect effects which may arise due to their finite length I. 

In the present contribution we investigate systematically the temperature depen- 
dence of the critical adsorption on a single spherical or rodlike particle, i.e., the case 
(t ^ 0, h = 0). (The generalization to the case h 7^ is straightforward but tedious.) In 
order to be able to treat spheres and cylinders in a unified way within a field-theoretical 
approach and for general spatial dimensions D it is helpful to consider the particle shape 
of a generalized cylinder K [30] with an infinitely extended 'axis' of dimension 5. The 
'axis' can be the axis of an ordinary infinitely elongated cylinder (5 = 1), or the midplane 
of a slab (5 — D — 1), or the center of a sphere (5 = 0). For general integer D and 5 the 
explicit form of K is 

K = jr = (r ± ,r|| ) e K^ 5 x M 5 ; |rj < (1.1) 

with r_|_ and rii perpendicular and parallel to the axis, respectively (see Fig.[l]). Note that 
r_|_ is a d- dimensional vector with 

d = D-5. (1.2) 



The radius R of the generalized cylinder K is the radius in the cases of an ordinary cylinder 
or a sphere and it is half of the thickness in the case of a slab. For the slab the geometry 
reduces to the much studied case of (two decoupled) half spaces. The generalization 
of D to values different from three is introduced for technical reasons because D uc = 4 
marks the upper critical dimension for the relevance of fluctuations of the order parameter 
leading to a behavior different from that obtained from mean-field theory valid for D = 4. 

In Sec. we discuss the general scaling properties of the local order parameter pro- 
files for critical adsorption on spheres and cylinders, in particular the behavior close to 
the particle surfaces and for small particle radii, respectively. In Sec. |T1] we consider the 
corresponding properties of the excess adsorption. In Sec. [TV] we present explicit results 
both for the order parameter profiles and for the excess adsorption in mean-field approx- 
imation. Section contains our conclusions. In Appendix [X| we discuss the two-point 
correlation function near a microscopically thin 'needle' at criticality. In Appendix [B| we 
determine a universal amplitude and a universal scaling function as needed in Secs.|IJ and 
|Tj , respectively. In Appendix 0, finally, we consider the general curvature dependence of 
the excess adsorption. 
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subspace || (dimension 8 = D-d) 




FIG. 1. (a) A generalized cylinder K with d < D (i.e., 6 = D — d > 0, see Eq. (|T7T|) ) and (b) a 
sphere as the special case of K with d = D as examples of particles with curved boundaries. D 
is the spatial dimension. The point r = (rj_, ry ) at which, e.g., the order parameter is monitored 
is also shown. 



II. ORDER PARAMETER PROFILES 

A. General scaling properties 

Asymptotically, i.e., in case of an infinitely large surface field, critical adsorption on 
the surface of a sphere or cylinder with radius R is characterized by an order parameter 
profile ($(r)) t which close to T c and for h = takes the scaling form 

(Hr)) t = a\tfP ± (s/Z ± ,R/Z ± ) (2.1) 

for radial distances s = r± — R > o from the surface larger than a typical microscopic 
length u. Here ( )t denotes the thermal average in the presence of a sphere or a cylinder 
and (3 is the standard universal bulk critical exponent. The scaling functions P± depend 
on two scaling variables: 

x± = s/£± , y± = R/i± • (2.2) 
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The scaling functions P± are universal while the nonuniversal bulk amplitudes a and £ 
are determined by the value = a|£| f° r £ < 0, t — > CT of the order parameter in 

the unbounded bulk and by the amplitudes of the correlation length, respectively. (In 
this work we are always concerned with the Ising universality class.) We note that the 
form of P± depends on the definition used for the correlation length £±. For definiteness 
we adopt £± as being the true correlation length fixed by the exponential decay of the 
two-point correlation function in real space pi]] . Bearing this in mind one finds that 
P +)b = P + (oo,y + ) = 0, P_ 6 = P_(oo,y_) = 1, and P±(x± -> oo,y±) - P± )b ~ exp(-ar±) 
where the prefactor in front of the exponential may by an algebraic function of x± (see, 
c.f., Sec.gVp. 

On the other hand, in the limit y± — > oo with x± fixed the scaling functions P±(x±,y±) 
reduce to 

P±( x ±) = P±{ x ±i°°) ) half-space , (2.3) 

corresponding to the half-space bounded by a planar surface. In recent years the scaling 
functions P±(x±) for the half-space have been intensively studied theoretically PJI^-IlI 



and compared with experiments P,0-|IB. For later reference we quote some of their 



properties. First, we note the asymptotic behavior ||T2 

P±(x±^0) - c±xZ P/v (2.4) 
x 1 + a± x± u + a' ± x± + b± x± + . . . 

with universal amplitudes c±, a±, a'±, and 6-t [^] which depend, however, on the defini- 



tion of the correlation length ||31|| . The ellipses stand for contributions which vanish more 
rapidly than x±. 



FIG. 2. (a) Universal scaling function P + (x+) for critical adsorption on a planar surface 
in spatial dimensions D = 4, 3, and 2. The estimates for P + (x+) in D = 3 labelled 'MC 
and 'RG' are obtained by Monte Carlo simulations and field-theoretical renormalization-group 
calculations. We display the corresponding data presented in Ref. [15]. The inset shows a scaled 
version of P+(x+) which reflects, in particular, the asymptotic behavior of P+(x+ — > 0). The 
open symbols indicate the corresponding values of c+ (see Eq. ( |2.4| ) and Table ||). (b) Same 
representation as in (a) for the universal scaling function P_(x_). The inset shows a scaled 



version of P_(x_) which contains the function U(x—) defined in Eq. (2.10) in order to facilitate 
a similar representation as in the inset of (a). 
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The exponent (3/u is the bulk 'scaling dimension' |33[ of the order parameter $. For 
convenience of the reader we quote their numerical values [53,3] : 



P(D = 4) = 1/2 , f3(3) ~ 0.328 , f3{2) = 1/8 , (2.5a) 

u{D = 4) = 1/2 , i/(3) ~ 0.632 , i/(2) = 1 , (2.5b) 

so that 

{(3/v) (D = 4) = 1 , (/?/!/) (3) ~ 0.519 , (2) = 1/8 . (2.5c) 



The leading power law in Eq. (|2.4j ) is a consequence of the fact that due to the presence 
of the symmetry breaking surface field the limit t — > must lead to a nonvanishing order 
parameter profile at t = so that 

($(r)) hM = = ac ± {s/it)-P /v ■ (2-6) 

(The subscript 'hs' stands for half-space). Thus the amplitude of the power law s~^/ u in 
Eq. ( [2.6D consists of the combination a (Co )^ °f nonuniversal bulk amplitudes and of the 
universal surface amplitude c±. Because the lhs of Eq. (|2.6|) does not depend on how the 
limit t = is approached, the universal amplitudes c + and c_ are related via the universal 
bulk ratio foVfo : 

The terms proportional to a± and a' ± in Eq. ( p.4|) correspond to regular contributions for 
t — > of the order parameter since 

($(r)) hM -> ($(r)) hM = [l + At + A'^ + ^iltl 2 - + ...] , t->0, (2.8) 
where A = ±o i ±(s/^q) 1 ^ 1 ' and A' = a ± (s/^ fc ) 2 / 1 ' are independent of the sign of t [12]. The 



term proportional to 6-t in Eq. ( |2.4|) corresponds to the first singular contribution p4| for 
t — > of the order parameter leading to the term B±\t\ 2 ~ a in Eq. (|2.8| ) with the bulk 
exponent a = 2 — z/L> and B± = b±(s / £,q) d ■ 

Figure summarizes theoretical results for P±(x±) for the spatial dimensions D = 4, 
3, and 2. In D = 4 the functions P± are given by the following mean-field expressions: 
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The results for D = 3 represent recent Monte Carlo (MC) simulations [|13| , |i"5f and field- 
theoretical renormalization group (RG) calculations fT2] , |r5|j . For D = 2 exact results are 
available from the semi-infinite two-dimensional Ising model ||35|| . Note that P±(x±) for 



any fixed value of x± decreases with decreasing D | 3q| . This reflects the general trend 
that critical fluctuations, which reduce the mean value of the order parameter, are more 
pronounced in lower spatial dimensions. 

The curves in the inset of Fig. ^ are scaled so that they demonstrate the leading 
behavior of P±(x± — > 0). (A similar representation will be used in Sec. IV.) In Table [j] we 



quote the corresponding numerical values of the surface amplitudes c± and a± according to 
Eq. ( |2.4| ). In order to achieve a presentation in the inset of Fig. 0(b) which reflects both the 
leading behavior of P_(x_ — > 0) and the exponential decay P_(x_ — > oo) — 1 ~ exp(— xJ) 
we introduce the function 



,2 



U(x-) = tanh — . (2.10) 

\X- + 1 / 

The curves in the inset of Fig. 0(b) comply with the above condition since U(x_ —>■ 0) ~ 
x 2 and t/(x_ — > oo) — 1 ~ exp(— 2xJ) so that in both limits the leading asymptotic 
behaviors of P_(x_) are not changed. 



TABLE I. Numerical values of the universal amplitudes c± and a± (see Eq. (2.4)) 



D 



C- 



4 

3, MC a 
3, RG' 3 

3, interpolation 
2 



v 7 ^ 1.414 

0.866 

0.717 
0.94 ±0.05 

0.803 



2 
1.22 
1.113 
1.24 ±0.05 
0.876 



-0.167 



-0.389 



-0.5 



i ~ 0.083 



0.129 



0.25 



a Reference [15] 

b Reference [^5| ; for a± we use in addition Eq. (48) in Ref. |jl2| for e = 1 

c Reference ]l^|; obtained from interpolating between the results in D = 4 — e and D = 2 
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In the remaining part of this section we discuss the new features of P±(x±,y±) which 
arise for y± < oo. First, we consider the short distance behavior for s — > 0, which 
corresponds to the limit x± — > with y± fixed. Then, we consider cases in which the 
radius R is small compared with £± as well as the distance between the sphere or the 
cylinder and the point for which the order parameter is monitored. This corresponds to 
the limit y± — > with x± fixed. It turns out that in this limit the behaviors for a sphere 
and a cylinder in D = 3 are qualitatively different. 

B. Short distance expansion 



The same reasoning leading to Eq. (|2.6| ) in case of the half-space yields for the gener- 
alized cylinder 

($(r)) t = = aCt( 7 ) {sKt)- p/v , 7 = s/R, (2.11) 

i.e., the universal amplitude c± in Eq. ( |2.6|) is generalized to the universal amplitude 
function C±(t) of the scaling variable 7. This function appears in the asymptotic behavior 

P±(x±,V± = x±/i) - C±( 7 ) xl p,v , x ± - 0, 7 fixed , (2.12) 

which underscores that C±( r y) is universal but depends on the definition of the correlation 
length []3T| (compare Eq. ( |2.4| )) and on the geometry. Since the limit 7 — ► must reproduce 
the behavior for the half-space one has C±(0) = c±. According to Eq. ( |2.7| ) the functions 
C+( 7 ) and C_( 7 ) are proportional to each other with 



C-{l) _ c_ (it 



>0 



P/v 

(2.13) 



Therefore it is sufficient to study only one of these functions, say, C + ( 7 ), which according 
to Eqs. ( gTjD and O can be written as 



C+( 7 ) = c + ($(r)) t = / (<K(r)) hM = . (2.14) 

In the case of a sphere, i.e., al = D, the universal scaling function C + ( 7 ) is known exactly 
for any spatial dimension D of interest by means of a finite conformal mapping from the 
half-space [[HJ]. It takes the simple analytic form 

C+( 7 ) = c+(l + ^j , rf = £>, (2.15) 
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and depends on D only via the corresponding values of c + and (3 /v. For large distances 
from the sphere, i.e. 7 ^> 1, the function £+(7) in Eq. (P-15|) decays as ^~@l v (compare 
the following subsection). In the opposite limit 7 — - > 0, which is of interest in the present 
subsection, one has 

^-°) =1 -^ T+ W/^ +0( ^), i=D . (2.16) 

We note, first, that an expansion such as in Eq. (|2.16| ) is expected to hold not only 
near the surface of a sphere but also near a smoothly curved surface of more general 
shape. According to differential geometry up to second order in curvature the short 
distance expansion near a (D — l)-dimensional surface of arbitrary shape involves only 



the geometric invariants K m , K£, and Kq with p7| , p8 

D-l 



K 

m 2 ^ R 

i=l 

and 

Kg = £ ^- (2.17b) 

pairs J 
i<j 

where Ri are the D — 1 principal local radii of curvature. In D = 3 Eq. ( [2.17|) yields 
the familiar expressions for the local mean curvature and the local Gaussian curvature. 
Second, we note that Eq. ( |2.16| ) reflects a property of the fluctuating order parameter field 
$(r) (or 'operator') in the outer space of the sphere itself. In this spirit also near a surface 
of more general shape the short distance expansion can be formulated in 'operator form': 

$(r)/($(r)) hM = (2.18) 
= [l]{l + K\ K m s + K 2 Kls 2 + k g K g s 2 + 0(s 3 )} + 0(s D ) . 

Here s is the distance of r to the nearest point of the surface. The terms on the rhs 
of Eq. (|2.18 ) should be interpreted as operators which are located at this point of the 



surface, i.e., as surface operators. On the rhs of Eq. (2.18) only such terms are shown 



explicitly which are proportional to the unity operator [1]. The corrections emerge from 
curvatures of higher order and surface operators of different types which are expected 
to scale with powers of s at least of the order s D (compare Eq. (|2.4|) for the half-space 
and Ref. pi|). Thus for D > 2 the terms displayed in Eq. ( |2.18| ) represent the leading 
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contributions of the short distance expansion of <3>(r). The dimensionless coefficients K\, 
k 2 , and kg depend on D but not on the shape of the boundary surface. Comparison of 
Eq. (p.18 ) with Eq. ( 12.16D for the sphere, or direct calculation near a surface of arbitrary 
shape flOf , yields 



D 



;2.i9) 



In D = 3 the coefficients k 2 and kq in Eq. ( |2.18| ) cannot be determined by comparison with 
Eq. fl2.16|) for the sphere alone. To this end one would need, in addition, the knowledge 
of the expansion for at least one curved surface of different shape (e.g., the surface of a 
cylinder). However, since in D = 3 the sphere has the property Kq = the comparison 
of Eq. ( |2.18| ) with Eq. ( |2.16| ) determines at least the sum 



K-2 + K G 
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D = 3 . 



(2.20) 



In Sec. [TV] we confirm Eq. ( [2.18 ) and determine k 2 , and k g for D = 4. 

The expansion ( [2.18| ) holds upon inserting it into thermal averages if the distance s 
to the curved surface — albeit being large on the microscopic scale — is much smaller 
then other characteristic length scales such as the correlation length or the distances to 
the remaining operators in correlation functions. For certain thermal averages such as 
the profile ($(r))t additional regular terms can occur on the rhs (compare Eq. ( |2.8| ) for 
the half-space). For spheres and cylinders, in particular, the expansion ( |2.18| ) determines 
the leading contributions to the scaling functions P±(x±,y±) generated by the surface 
curvature. For the surface of a generalized cylinder the curvatures in Eq. (|2.17|) are given 
by 

d-l 1 
2 R 



(2.21a) 



and 



Kc 



(d-l){d-2) 1 



2 R 2 ' 

By inserting these expressions into Eq. (|2.18|) and using Eq. ( [2.1|) one obtains 



(2.21b) 



P ± (x ± -> 0,y± 



+ 



«2 



sdc 



C± X± 



1 + Kl 



d-l 



■7 



'2.22) 



(d-l){d-2) 



7 2 + 0(7 3 
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for x± — > with y± fixed so that 7 = x±/y± — > (the subscript 'sde' refers to the short 
distance operator expansion in Eq. ( j2.1S| )). 

Note that Eq. (|2.22| ) does not contain the above mentioned regular contributions. Upon 
employing Eq. (|2.4j ) for the half-space, however, one can obtain the asymptotic expansion 
of P±{x± — > 0, y-t) including the leading regular term. Assuming that for t ^ the profile 
($(r)) i is still analytic in 1/R one finds 



P±(x±^0,y±) 



-P/v 



+ 



(d-iy 

«2 T h 



(d- l){d-2 



d-1 
1 + «i — n — 7 



(2.23) 



2 1 - i/f 1 
7 + a± x ± + . . . 



for x-t — > with y± fixed so that 7 = x±/y± — ► 0. The ellipses stand for contributions 
which vanish more rapidly than 2^ if D > 2 []39| . The universal amplitudes c± and a± are 
the same as in Eq. ( |2.4p (see Table |) and «i, K2, and kg are from Eqs. fl2.19| ) and ( |2.20 ). 
In Sec. |V| we shall confirm Eq. (|2.23|) explicitly for D = 4. 



C. Spheres and cylinders with small radii 

In this subsection we consider spheres and cylinders whose radii are much smaller 
than other characteristic lengths such as the correlation length or the distance between 
the particle and the point at which the order parameter is monitored. In these limiting 
cases the effect of the particle upon the fluctuating order parameter distribution can be 
represented by a 5-function potential located at the center of the particle which enhances 
the value of the order parameter. It is instructive to consider this expansion for the 
generalized cylinder K for which this 5-function potential is smeared out over its axis, i.e., 
the Boltzmann weight exp(— STix) of K, where STix is the difference of the Hamiltonian 
describing the system with and without the presence of K (whose axis includes the origin), 
can be systematically expanded in a series with increasing powers of R p!9|j23[], i.e., 



exp{-5H K ) oc 1 + S d>D R^ v - D+d w K + ... (2.24) 



where £^d is an amplitude and 



/ d 8 r\\ $(r_L = 0,r|j) , d < D , 
w K = { RS (2.25) 
$(0), d = D. 



13 



Here only the leading nontrivial contribution for i? — f is shown explicitly and the ellipses 
stand for contributions which vanish more rapidly for R — > 0. For the case that K is a 
sphere i.e., d = D, the amplitude £d,d is equal to the ratio A® / , where A® and 
B$ are amplitudes of the half-space profile 

($(r)) hM = = A*(2s)-^ (2.26) 

at the critical point of the fluid for the boundary condition j corresponding to the critical 
adsorption fixed point and of the bulk two-point correlation function 

($(r)$(0)) M = = B^r' 2 ^ (2.27) 



at criticality, respectively. The ratio (A®) / B$ is universal |23|. The comparison of 
Eq. ( |2.26| ) with Eq. ( |2.6|) yields the relation 

A* = a(2Ctf /v c ± (2.28) 

between the nonuniversal amplitudes A%, a, ^ and the universal amplitude c±. 

It is crucial to observe that Eq. ( |2.24! ) is only valid if the exponent of R is positive, i.e., 

P/v - D + d > . (2.29) 

Figure shows as a dashed line d = D — ((3/u)(D) ||(| in the (d, _D)-plane which sepa- 
rates generalized cylinders K which are relevant perturbations for the fluctuating order 
parameter field (such as the strip in D = 2 or the plate in D = 3) from those which 
are irrelevant and for which Eq. (2.24) applies. Thus, the disc in D = 2 and the sphere 



in D = 3 represent irrelevant perturbations whereas the cylinder in D = 3 represents a 
relevant perturbation. This implies that an infinitely elongated cylinder in D = 3 gen- 
erates a perturbation of the order parameter from its bulk value whose spatial extension 
is only limited by the bulk correlation length. The order parameter profile becomes even 
independent of R in the formal limit R — > 0, i.e., if the cylinder radius R becomes mi- 
croscopically small ^TJ. The critical adsorption transition on such a microscopically thin 



'needle' is characterized by critical exponents which need not be equal to the correspond- 
ing exponents for the bulk or the half-space (see Appendix [A]). For the disc in D = 2 and 
the sphere in D = 3, in contrast, the deviation of the order parameter from its bulk value 
vanishes in the limit R — > 0. This is reflected by Eq. Q2.24 ) which is — apart from the 



condition in Eq. (|2.29 ) — only valid if the small radius R is still large on the microscopic 
scale. 
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FIG. 3. Diagram of generalized cylinders K which behave — in the renormalization group 
sense — as relevant or irrelevant perturbations of a fluid near criticality. The parameter 
d < D characterizes the shape of K and D is the space dimension (see Eq. The point 

(d, D) = (2, 2) corresponds to a disc in D = 2 and the points (3, 3) and (2, 3) to a sphere and 
an infinitely elongated cylinder in D = 3, respectively. The line with D = 4 and arbitrary d 
represents the upper critical dimension for which the mean-field results for the adsorption pro- 
files are exact (see Sec.|TV|). The open circles indicate points (d, D) for which d = D — ((3/v)(D) 
within the Ising universality class. These points are connected by the dashed line so that within 
the shaded region above it the small radius expansion fl2~24l) is valid (see Eq. (]2~29|) ) and K 
represents an irrelevant perturbation. Points (d, D) below the broken line, such as the cylinder 
in D = 3, are characterized by the fact that the order parameter at large distances from K 
deviates from its bulk value even if the radius R is microscopically small, which means that K 
represents a relevant perturbation. 



The line d = D — {(3/v){D) itself corresponds to marginal perturbations leading to a 
behavior which in general is different from Eq. ( |2.24| ). We shall neither discuss this nor 
the crossover from marginal behavior to the behavior described by Eq. fl2.24p which may 
arise for points closely above the line. The line d = D — (f3/v)(D) includes, in particular, 
the generalized cylinder with d = 3 and D = 4 (see Sec. fV\). The marginal behavior 
of this particular generalized cylinder, however, is not typical for spheres or cylinders in 
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D = 3, which are strictly irrelevant or relevant perturbations, respectively. Therefore in 
the following we shall generally speak of spheres if d = D and speak of cylinders if d = 2 
and D > 3 (see Fig.|). 

We now turn to the consequences of Eq. ( ]2.24j) and the related properties of the uni- 



versal scaling functions C+(j) and P±(x±,y±) defined in Eqs. ( P-14|) and (|2.1|) . First 



we consider the sphere, i.e., d = D, for which Eq. ( j2.24| ) applies. In order to obtain the 
asymptotic behavior of C + (j) for 7 = s/R — > 00 one basically replaces ($(r)) 4 = in the nu- 
merator on the rhs of Eq. (|2.14j) by the bulk two-point correlation function ($(r)$(0))&.f = o 
where r can be replaced by s in leading order. By using Eqs. ( |2.26| ) and ( |2.27|) in con- 
junction with £d,d = A®/B<$, one finds 

C + (j^oo) - c + (j) , d = D. (2.30) 

This checks with the exact result for £+(7) in Eq. (|2.15|) . For t ^ and y± = R/£± — > 
with x± = s/£± > fixed the same steps as above lead to 

^±(^±-0) - P ±tb ^ c ± (2y ± )^xZ 2l3/u F±(x ± ) , d = D . (2.31) 

Here T± is the bulk universal scaling function defined by (compare Appendix ||) 

($(r)$(0)>£ = B„r~ 2 ^ ^±(r/e±) (2.32) 

where the superscript C denotes the cumulant of the correlation function. Equations 
( j2.32j ) and (|2.27D imply the normalization ^±(0) = 1. According to the definition 
of £± as the true correlation length (compare the discussion below Eq. ( |2.2D ) one has 
J r ±(x± — > 00) ~ exp(— x±). Equation (|2.31|) implies that P±(x±,y± — > 0) decays to its 
bulk value with the power law ~ y±" for any value of x±. 

Next, we consider the cylinder, i.e., d = 2 and D > 3. In accordance with the 
discussion above this object represents a relevant perturbation. The order parameter 
deviates from its bulk value even for R — > so that the universal scaling functions £+(7) 
and P±(x±, y±) remain finite in the limit 7 = s/R — > 00 and y± = R/£± — > 0, respectively. 
For t — this implies 

C±(7 -> 00) -> n± , d = 2, D > 3 , (2.33) 

which defines new universal amplitudes n± with n_/n+ = (£q" '/ £0)^^ ■ Equations (|2.33|) 
and ( |2.11| ) imply that for t — the order parameter profile near a thin 'needle' at distances 



s from the needle large compared with microscopic lengths takes the form 
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($(r)) t=0 = an + (s/ti)-^ , 'needle' . 



(2.34) 



For t ^ one finds 



P±(z±,y±^0) - P±(x±,0) = N ± (x ± ) , rf = 2, D>3 



(2.35) 



where the new universal scaling functions N±(x±) characterize the critical adsorption 
profile on a thin needle. In the limit x± — > they behave as 

N±(x± -»• 0) -> -n±xg P,v , 'needle' , (2.36) 



with n± from Eq. (|2.33|) (compare Eq. (|2.4j) for the half-space). In Sec. IV we shall confirm 
Eqs. (|2.33| ) - ( |2.36| ) and calculate n± and N±(x±) explicitly for D — 4. Figure [| summarizes 
the various types of limiting behavior of the scaling functions P±(x±,y±). 



(2.6) 

\ (2.4) (2.3), half-space 

CO 







(2.23) 



(2.16) 



(4.14) 




(2.12), slopes 1/7 



30), (2.33) 



(2.31), (2.35) 



bulk 



(2-36) 



oo 



FIG. 4. The universal scaling functions P±(x± = s/£±,y± = R/£±) exhibit distinct behaviors 
for various limits of the variables x± and y±. The position of an equation number in the 
(x, y)-plane indicates the corresponding limiting case to which it applies. The case y± = oo 
corresponds to the half-space (see Eq. fl2.4|) and Table |). The short distance behavior for 
s — > corresponds to the limit x± — * with y± fixed (see Eq. ( 2.23[ )). The universal functions 
C±(7 = s/R) in Eq. ( |2.12| ) characterize the behavior for £± — > oo, i.e., the limit x± — * with 
fixed 'slope' 7 _1 = y±/x±. The behavior for small i? corresponds to the limit y± — > with x-t 
fixed (see Eqs. (jOHD and ( p5| )) 
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III. EXCESS ADSORPTION 



A. General scaling properties 



Close to T c the total enrichment of the preferred component of the fluid near the 
surface of the generalized cylinder is proportional to the excess adsorption T(t, R) defined 



as 



r(t,R) 



d D r 



f, b 



(3.1) 



v 



I s fl d I ds (s 



R) d - L ($(r)) t - ($) 



t,b 



Here V = M. D \ K is the volume accessible to the critical fluid, i.e., the total space except 
for the volume occupied by the sphere or the cylinder K, and = 2n d ^ 2 / T(d/2) is the 
surface area of the d-dimensional unit sphere. For a sphere one has 5 = D — d = and 
I s = 1 whereas for a cylinder d = 2 so that in D = 3 the quantity I s = I is the length of 
the cylinder. In order to obtain the scaling behavior of T(t, R) we split the s integration 
into the intervals < s < a and s > a where o represents a typical microscopic length 
| 14j| . For the second interval one can use Eq. ( |2.1D which leads to 



r{t,R) 



I s n n 



I(t;a,R) + J±(y±;<r/Z±) 



with 



I(t;a,R) = J dsis + R)*" 1 ($(r)) t - 
o 



(3.2) 



(3.3) 



^±(y±;o-/C±) = a\tf£± J dx±(x± + y±. 

<r/t± 



d-l 



P±{x±,y±) - P±,i 



(3.4) 



In the limit t — > the integral in Eq. ( |3.3| ) remains finite and yields a nonuniversal constant 
which is subdominant to the diverging contribution J± in Eq. ( |3.4| ). In order to clarify 
the dependence of the latter on y± = R/£± we decompose it according to 



(3.5) 



x < / dx± 



P±{x±)-P±, 



18 



where P±(x±) = P±{x±, oo) are the scaling functions for the half-space and 



(3.6a) 



x± 

y± 



+ 1 



d-l 



P±{x±,y±)-P± >b - P±{x ± )-P ± , 



The comparison of Eq. Q3.5| ) with Eq. ( |3.2p shows that the first term in curly brackets in 
Eq. ( |3.5|) renders a contribution of the interval s > a to the excess adsorption per unit 
if the surface of the particle would be planar times the area A = l s VtdR d 1 of 
the actually curved surface of a sphere or a cylinder. Thus the functions G± in Eq. (|3.5|) 
reflect the deviation of the excess adsorption on a curved surface from that on a planar 
surface beyond pure geometry. 

In order to reveal the behavior of J±(y±; & / £±) for t — > it is necessary to study 
the integrals in Eqs. ( |3.5|) and ( |3.6a|) at their lower bounds a/£± — > 0. The integral in 
Eq. ( |3.5| ) can be analyzed along the lines of Sec. IIA in Ref . . One finds that for 
D < 4 its contribution to T(t —>■ 0, R) leads to the power law singularity ~ \^f~ v I [y — 0) 
corresponding to a planar surface. The proper limit D 4, for which v — > f3 (see 
Eq. ( |2.5|) ), is accomplished by the presence of a term constant with respect to t which 
also diverges for D /* 4 such that the sum leads to a contribution to T(t — > 0, R) which 
diverges logarithmicly in D = 4. The integral in Eq. ( |3.6aj ), however, remains finite for 
a/£± — > also in D = 4 since the singular behavior for x± — * of the first term in curly 
brackets is cancelled by the second term. This implies that the function 



G±(y, 



G±(y±; <r/i± = 0) 



(3.6b) 



does indeed represent the leading behavior for t — > 0. Note that G±(y±) is universal 
because it depends only on the universal scaling functions P±(x±,y±) and P±(x±). In 
sum one finds 



T(t^0,R) - Aog{ fa |f| ^ 1 + \tf- v G ± (y±)} 
with the universal numbers (see Table ||) 



D< 4 



9± 



(i/-/3) [p±(x ± )-p ± 



V c± 



D < 4 , 
D = 4, 



(3.7) 



(3. 
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where the second line is the limit for D / 4 of the first line. According to the above 
discussion one has Cr±(oo) = 0. Equation fl3.7| ) generalizes the corresponding Eq. (2.10) 
in Ref. [14[] for a planar surface by the additional second term in Eq. fl3.7|) . A quantity 
accessible to experiments is the ratio 

r(+\t\,R) 



K*(\t\,R) 



\t\,R) 



(3.9a) 



of the excess adsorptions above and below the critical point. The leading behavior of IZq, 
for \t\ — > is characterized by a universal function R$(\t\, y+, y~), i.e., 



K${\t\^0,R) -> R*(\t\,y+,y-) 



(3.9b) 



which can be read off from Eqs. ( |3.9a|) and 
to the universal number [|T4| 



For a planar surface this function reduces 



Ri 



— — , halt-space . 
£o 9- 



(3.10) 



Table U summarizes theoretical results for g± and _R<i> corresponding to the half-space. For 
the curved surface of a sphere or cylinder, however, we shall see that for y± = R/£± — > 
the divergence of the second term in curly brackets in Eq. ( |3.7|) is more pronounced than 
the divergence of the first term. In order to clarify this aspect we now discuss the behavior 
of G±(y±) in the limit for large and small values of y±, respectively. 



TABLE II. Numerical values of the universal numbers g± and i?$ (see Eqs. ( |3.8j ) and ( |3.10| )). 



D 



9+ 



9- 



4 

3, MC a 
3, RG a 

3, interpolation 13 
2 



V2/2 ~ 0.707 
0.663 
0.581 
0.69 ±0.1 
0.910 



1 

0.599 
0.438 
0.56 ±0.1 
0.0818 



2.28 ±0.1 
22.236 



a Reference [15] 

b Reference [14|; obtained from interpolating between the results in D = 4 — e and D = 2 
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B. The scaling function G±(y±) for y± — > oo 

For i/i > 1 we assume that G±(y±) is analytic in y± x = £±/R so that it can be 
expanded into a Taylor series around y± x = 0, i.e., 

G ± (y ± - oo) = a^^" 1 + 6±^ 2 + . . . (3.11) 

(recall Cr±(oo) = by definition). The coefficients a^ D and 6^ are dimensionless and 
universal and depend on the space dimension D and on the shape of the generalized 



cylinder, i.e., on d. The validity of the expansion (|3.11|) is plausible since in the limit 



£±/R <C 1 the thickness ~ £± of the adsorbed layer is much smaller than the particle 
radius R so that a small curvature expansion should be applicable to T(t — > 0,R) in 
Eq. ( |3.7| ) in which a surface term ~ I s R d_1 is followed by successive terms ~ R d ~ 2 , R d ~ 3 , 
etc., generated by the surface curvature. The terms on the rhs of Eq. ( |3.11| ) correspond 
to the leading curvature contributions to this expansion and imply that for decreasing 
values of y± = R/£± ~ \t\ u the leading corrections to the first term in curly brackets in 
Eq. ( |3.7| ) — corresponding to a flat surface — are of the order \t\^~ u \t\~ ni/ with n = 1,2, 
etc., and start to dominate it as soon as y± < 1. 

Similar as for the short distance expansion of the order parameter (compare Eq. ( |2.18| )) 
the first Taylor coefficients of the expansion of G±(y± — ► oo) also determine the curvature 
parameters of a particle K, of more general shape provided its surface S is smooth and all 
principal radii of curvature are much larger than £±. In this case one expects an expansion 
of the general form 

T(t^0,R) = (3.12) 

m^-i 



X 

s 



where the curvatures K m and Kq are given for a (D — l)-dimensional surface of general 
shape in Eq. (|2.17 ). For the special case that the particle K, is a generalized cylinder K 



the curvatures K m and Kq are given by Eq. ( |2.21| ) and the comparison of Eq. ( |3.12j ) with 
Eqs. O and (|37TT|) yields 



= W- V a± e ± (3.13a) 



and the relation 
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{d-lf , x± (d-l)(d-2) 



+ A 



G 



(3.13b) 



The curvature parameters Aj 4 , A^, and A^ depend on D but should not depend on the 



shape of K, i.e., on d in Eq. ( |3.13|) . This implies a corresponding dependence on d of 



,D 



and 6^£), which provides an important consistency check for the validity of Eq. (|3.12|) . In 



Sec. [IV| we confirm this dependence and explicitly calculate the curvature parameters for 
D = 4. 



C. The scaling function G±(y±) for y± — ► 



In order to investigate this limit it is convenient to consider the function 



(3.14) 



dx ± \{x ± + y ± ) d - 1 P±(x ± ,y ± )-P ± 



d-l 

y± 



P±(x±) - P±, 



Because for y± — > the function Pi^iiZ/i) m Eq. Q3.14D behaves qualitatively different 
for spheres as compared with cylinders (see Sec. [II C|) we treat them separately. 

First, we consider the sphere, i.e., d = D, for which Eq. ( [2.31|) holds for the behavior 
of P±(x±,y± — > 0) if x± > is fixed. However, the integral in Eq. ( |3.14| ) starts at the 
lower bound x± = where this condition for x± is violated. Therefore we temporarily 
split the x± integration into the intervals < x± < ^fy± and x± > ■ s /y±. In the latter 
interval the relation s/R > y± 1 ^ 2 — > oo holds so that Eq. ( p. 31] ) is applicable. This leads 
to (the curly brackets correspond to those in Eq. ( |3.14| )) 



/y± 



y±G±(y±) 

oo 

+ J dx± 



dx4 



(3.15) 



,x± + y±, 



>D-1 



c± (2y ± )^ 'x ± W/l 'F ± {x ± 



yf- 1 



P±(x±)-P±, 



In the second integral the variable y± in {x± + y±) D ~ l can be replaced by zero in leading 
order for y± — > and the term proportional to y±~ l , i.e., the half-space contribution, can 
be dropped since D — 1 > (3/ v. The resulting integrand is integrable for x± — > so that 
the lower bound Jy± can be replaced by zero and one can readily show that the resulting 
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second integral in Eq. (|3.15 ) dominates the first integral for y± — > 0. This leads to the 
final result 

G±(y± - 0) - u± y± D+1+ ^ , d = D, (3.16) 
with the universal amplitudes 

oo 

u ± = c±2^ u J dx±x±- l - W,u F±{x±) . (3.17) 



The numerical values of u± for several spatial dimensions D as presently available are 
summarized in Table |TTT] (see Sec. |IV B| and Appendix ||). Equation ( |3.16| ) in conjunction 



with Eq. ( |3.7| ) renders the power law divergence 

I\t ->0,fl) ~ \t\ 2(3 - Du , d = D. (3.18) 



For t — > the universal function y + , y_) in Eq. (|3T9|) approaches the universal 

number 



i^(0,0,0) = lim 1) ' )7"( = ( 7zr ) — , d = I>. (3.19) 



t- 



r(+ 




,R) 


r(- 


t 


,R) 



Whereas the dependence of the function R$(\t\,y + ,y-) on y + , y_ is universal but de- 
pends on the definition used for the correlation length, we note that the limit _R$(0, 0, 0) 
in Eq. ( |3.19|) is independent of the definition for the correlation length. The reason is that 



according to the first parts of Eqs. ( |3.19| ) and ( |3.1| ) R$(0, 0, 0) can be expressed in terms 



of the order parameter profile ($(r)) t without resorting to the notion of the correlation 
length at all. (The same holds for the number in Eq. ( |3.10| ) JHJ].) 



TABLE III. Numerical values of the universal amplitudes uj± and v± (see Eqs. ( [3.17 ) and 
( 3.21] )) as presently available. 



D 


u + 




v+ 


V- 


4 


4y/2 ~ 5.657 


8 


1.90 


1.86 


3 


1.53 ±0.05 


~ 1.47 






2 


0.515 


0.0501 


not defined 


not defined 
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Next, we consider the cylinder, i.e., d = 2 and D > 3, which represents a relevant per- 
turbation (see Sec. [11 (J[ ). In this case Eq. ( |2.35| ) holds for the behavior of P±(x±, y± — > 0). 
Thus y± _1 G±(y±) = y±G±(y±) becomes independent of y± in the limit y± — > and tends 
to the constant given by the integral in Eq. (|3.14j ) with y± replaced by zero. This leads to 

G ± (y ± -> 0) -> v ± y± x , rf = 2,D>3, (3.20) 



with the universal numbers 



v± = J dx± x± 





W±(s±)-P±, 



(3.21) 



Numerical values of v± are only known for D = 4 at present (see Table |IT| and Sec. |IVB]) 



Note that the numbers v± are not defined for D = 2 (compare Sec. [II C|) . The exponent 
in the power law in Eq. ( |3.20|) differs from that in Eq. ( p,16| ) — apart from the difference 
generated on purely geometric grounds — by j3/v due to the different small y± behavior 
of the critical adsorption profiles for spheres as compared with cylinders. Equation ( |3.20| ) 
in conjunction with Eq. (|3.7f) leads to 



T(t -»• 0,R) ~ \t 



f3-2v 



d = 2 ,D > 3 . 



(3.22) 



For t — > the universal function R$(\t\,y + , yJ) in Eq. (|3.9[) tends to the universal number 



i2*(0,0,0) = lim 



n+\t\,R) . . (ti\ 2 v + 



°r(-\t\,R) U 



— , d = 2 ,D>3 



v. 



(3.23) 



with v± given by Eq. ( |3.21| ). Again, the limit R<& (0,0,0) is independent of the definition 
for the correlation length (compare Eq. ( |3.19 )). 



IV. MEAN-FIELD THEORY 



The critical fluctuations of the fluid are described by the standard Hamiltonian [[7].Q 



nm = J dvfyv®) 



2 + 1$ 2 

2 



u 
24 



(4.1) 



for a scalar order parameter field $(r) supplemented by the boundary condition $ = +oo 
at the surface of the sphere or the cylinder corresponding to the critical adsorption fixed 
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point f[2" |. The position vector r G M 15 covers the volume V = IR D \ K accessible to the 
critical fluid. The parameter r is proportional to t = (T — T c )/T c and u is the $ 4 coupling 
constant. The thermal average ($(r)) corresponding to the Hamiltonian in Eq. Q4.1|) with 



the boundary condition for critical adsorption at the surface of K can be systematically 
expanded in terms of increasing powers of u, i.e., 

($(r)) = x -\m(s;R,r) + 0(u)] , s = r ± - R . (4.2) 
V u i J 

The leading contribution to this expansion corresponds to the mean-field result for the 
order parameter profile which becomes exact in the limit D /* 4. The profile m(s;R,r) 
is determined by minimizing 7i{$} which leads to the Ginzburg-Landau type equation of 
motion 

m"(s) + — — m'(s) —Tm(s) = m 3 (s) , (4.3) 
which is supplemented by the boundary conditions 

m{s ^0;R,t) -> — , (4.4a) 



T > , 

m(s — > oo; R, r) — > m± i f, = ^ (4.4b) 

|r| 1/2 , r < . 

The parameter r is related to t — (T — T c )/T c and to the correlation length £± by 

- 1 / 2 , r > , 

\r\- 1 ' 2 , r < . 





& = ^i*r 1/2 = < 1 , ,_,„ „ (4-5) 



and (compare the text following Eq. ( |2.2| 



a | t |V2 = /6 | r |l/2 ^ T<Q (46) 



Equations ( |4.3|) and ( [4.4|) uniquely determine the profile m(s;R,T), which is related 
to the universal scaling functions P±{x±,y±) in Eq. ( |2.l|) by 

P±(^ ±) y±) =|r|- 1 / 2 m(s;fl,r) , D = 4 . (4.7) 
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Within the present mean-field approach this scaling form holds because the lhs of Eq. ( |4.7| ) 
is dimensionless so that it depends only on the dimensionless variables s|t| 1//2 and Rlr] 1 ^ 2 , 
where \t\ 1 ^ 2 is related to £± by Eq. fl4.5|) . Similarly, the universal scaling function £+(7) 
in Eq. (|2.11| ) is given by 

C + (-f) = sm(s;R,r = 0) , D = 4 . (4.8) 



Equations ( |4.7| ) and ( |4.8|) lead to exact results for D / 4 and generalized cylinders K 
with arbitrary d in the interval 1 < d < D (compare Fig. [3D. Equations Q4.3Q and (4.4) can 
be solved numerically, e.g., by means of a shooting method [f£2]| . In some cases analytical 
solutions are available. In the remaining part of this section we present the corresponding 
explicit results for the problems discussed in Sees. [TT] and |ITT| . 



A. Order parameter profiles 



We start with the universal scaling functions P±{x±,y±) according to Eqs. ( |4.7|) and 
( |4.3| ). Figure |] shows their behavior as function of x± — s/£± for various values of the 
parameter y± = R/£± in the case of a sphere (i.e., d = D with D = 4, compare Sec. |ll (J| 
and Fig.^). For our presentation we choose a scaled form which reflects the behavior of 
P±{ x ± ~ > 0, 2/±) (compare the inset of Fig. § here (3/v = 1). Accordingly, the curves start 
at the mean-field values c+ = a/2 and c_ = 2 (see Eq. Q2.23Q and Table[I]) , respectively, 
and decay exponentially for x± — > 00. 



FIG. 5. (a) Scaling function for a sphere in mean-field approximation (i.e., d = D 

with D = 4) as function of x+ for several values of y+. The curves start at the value c + (compare 
the inset in Fig. [2]). The curve for y + = 00 corresponds to the half-space profile (see Eq. (|2.9[ )). 
The curves decrease with decreasing y + and vanish in the limit y+ — * 0. The slopes at x + = 



can be read off from Eq. ( 4.16| ) with s/R = x + /y + . (b) Same representation for P_(x_,?/_) 
The function U(x—,y—) introduced for convenience is defined in Eq. 
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1 2 3 4 5 6 

x + = s/^ + 




x_ = s/t 
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x_ = s/t 
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FIG. 6. (a) Scaling function P+(x+, y+) for a cylinder in mean-field approximation (i.e., d = 2 
and D = 4) as function of x+ for several values of y+ (compare Fig. ||). The curves decrease with 
decreasing y + . In the limit y + — > they tend to the curve corresponding to P+ (x + , 0) = A^ + (x + ) 
describing the critical adsorption profile on a thin needle. The latter curve starts at the value n+ 
(filled symbol), (b) Same representation for P_(x_,y_). The function U(x_) equals U(x-,oo) 
from Eq. 



In order to deal with the rapid decay of P_(x_ — > oo, y_) — 1 for small values of ?/_ in the 
case T < T c we replace the function U(xJ) in Eq. ( |2.1(J| ) by 



U{x-,yJ) = tanh ( - X + - V ^ 1 j , d = P> = 4 , (4.9) 

so that Z7 (x__, oo) = C/(x_). We emphasize again that the function U(x-, y_) is introduced 
only in order to facilitate an appropriate representation of P_(x_,yJ) which reflects the 
behavior for both small and large values of X- for all values of y_ (compare the discussion 
related to Fig. 0). The overall dependence of P±{x±,y±) on the parameter y± is in line 
with Eq. fl2.31|) and the related discussion. Starting from the half-space behavior for 
y± = oo the profiles P±(x±,y±) — P± t b decrease with decreasing y± for any fixed value of 
x± and vanish in the limit y± — > 0. 

This latter behavior for a sphere differs from the corresponding one for a cylinder 
(i.e., d = 2 and D = 4, compare Sec. |II C| and Fig.Q). Figure ^| shows the corresponding 



behavior of P±{x±,y±) for various values of y±. For the presentation of P_(s_, yS) in the 
case of the cylinder the function U{xJ) = U(x-, oo) is suitable for all values of y_. The 
overall dependence on y± is characterized by the fact that P±(x±,y±) does not vanish 
in the limit y± — > but rather tends to the finite limit function P±(x±,0) = N±(x±) 
corresponding to the critical adsorption profile on a thin needle. This is in line with 
Eq. ( |2.35|) and the related discussion. According to Eq. (|2.36| ) the curves for y± = start 



at the mean-field values of the universal amplitudes n± which are given by 

n+ = 1 , n_ = V2, D = 4 . (4.10) 



Next, we consider the universal scaling function £+.(7) in Eq. (|4.8|) corresponding to 
r = 0. Figure shows its behavior for the present case D — 4 and several values of d. The 
function C+(y) starts at C+(0) = c + (see the text following Eq. ( |2.12|) ) and for 7 — » 00 it 
vanishes in the case of the sphere (i.e., d = D with D = 4) while it tends to the finite 
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D = 4 



half-space (d = 1 ) .. 



cylinder (d = 2) 




FIG. 7. Scaling function £+(7) in mean-field approximation (i.e., D = 4) as function of 
7 = s/R for several values of d (compare Fig. |3j) . All curves start at the same value C+(0) = c+ 
corresponding to the half-space but for large 7 the behaviors are qualitatively different. For a 
sphere £+(7 —* 00) vanishes as ^~^l v with (3jv = 1 whereas for a cylinder it tends to the finite 
value n + . In the marginal case d = 3 it vanishes only logarithmically. 



number n + from Eq. ( |4.10|) in the case of the cylinder (i.e., d = 2). In the marginal 
case d = 3 it vanishes only logarithmically (compare Sec. |11 C| and Fig. 0). For r = 
analytical solutions for Eq. ( f4.3| ) are available in some special cases even for d > 1. For 
example, for r = and d = 4 the differential equation can be solved by relating it to 



the generalized Emden- Fowler equation [ 22 1, or by realizing that it represents a so-called 



'equidimensional' equation (see, e.g., section 1.4 in Ref. [0). The result is 

2V2R 



m(s; R, t = 0) 



s (s + 2R) 



d 



(4.11) 



which agrees with Eqs. ( |4.8|) and (|2.15|) . It is worthwhile to note that an analytical solution 
also exists for the particular case r = and d = 5/2 for which 
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m(s;R,r = 0) = — / ' 2 d = 5/2 , (4.12) 

V ; s + R - y/R (s + R) 

from which the corresponding result for £+(7) can be inferred from Eq. ( |4.8|) . Finally, for 
t = 0, R = 0, and arbitrary values of d < 3 one finds 



/3 

m(s;# = 0,r = 0) = , s = r ± , d < 3 , (4.13) 



so that C + (oo) = n + = y/3 — d > for d < 3. This is in line with Eq. ( |4.12|) as well 



as with Fig. |3] and the related discussion, according to which in D = 4 a generalized 
cylinder K with d < 3 represents a relevant perturbation. In particular, Eq. ( |4.13| ) leads 
to m(rj_) = l/r_i_ representing the mean-field profile near a thin needle (for which d = 2, 
compare Sec. |II C|) at criticality. In Appendix [A] we calculate the corresponding two-point 
correlation function in the presence of this profile. 

The asymptotic behavior of P±{x±,y±) for D = 4 in various limiting cases can be 
derived directly from the defining Eqs. ( f4.7|) and (|4.3|) even for those cases in which no 
full analytical solution is available . 

(i) x± = s/£±. ^> 1 with y± = R/£± fixed. Because in this limit m(s;R,r) — m± t b is 
exponentially small, in Eq. fl4.3| ) one can neglect powers of it larger than one. This leads 
to a linear differential equation for m(s]R,r) — m± j f } which can be solved in terms of 
modified Bessel functions K a (x± + y±) and I a (x± + y±) with index a = (d — 2)/2. 
Here only the decreasing function K a must be considered. Using the asymptotic behavior 
of K a for large arguments [4~5| and Eq. ( |4.7| ) one finds 



d-l 



P±(x±,y±) ~ P±,b -> A±(y±;d) {x± + y±) 2 exp(-x±) , x± > 1 , (4.14) 
where the amplitude- function remains undetermined by the present method 



and must be evaluated numerically in general. An important feature of Eq. ( 4.14j) is that 
for increasing values of x± the exponential decay is enhanced by the algebraic prefactor. 
Thus the decay is weaker if x± y± and it is weaker for a cylinder (d = 2) than for a 
sphere (d = D with D — 4). This is expected because the perturbation of the bulk fluid 
is strongest in case of the half-space (corresponding to x± -C y±), less for a cylinder, and 
even less for a sphere. 

(ii) x± — > with R and ^± fixed, i.e., s —> (compare Sec. IIB ). The leading behavior 
of m(s — > 0; R, r) can be obtained by inserting the ansatz 
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m(s;R,T) = h a + a%s + G(s 2 ) (4.15) 



s 



into Eq. (|4.3|) and fixing the coefficients ao and a± such that the prefactors of the most 
singular powers of s cancel. The resulting expression can be arranged so that it takes the 
form implied by Eqs. ( |2.23| ) and ( |4.7| ), i.e., 

m(s ^0;i i ,r) = f (4-16) 



+ 



5 (d- l) 2 1 (d— l)(d-2)-i s 2 1 _, 



L9 



i? 2 6 TS + 



Equation (|4.16|) confirms that the leading dependence on r is analytic with the same 

prefactor as for the half-space. In addition, one can read off the curvature parameters 

according to mean- field theory, i.e., 

15 1 
^ = -3, K 2 = -, k g = --, D = A. (4.17) 

We note that the derivation of Eq. (|4.16|) implicitly contains a consistency check for the 
validity of the general curvature expansion (|2.23|) . For example, the parameters K2 and kq 
are fixed by considering only two different types of generalized cylinders K with nonvan- 
ishing curvature, e.g., those for which d = 2 and d = 3. Equation ( 4.17] ), however, holds 



for any value of d in the interval 1 < d < D with D = 4 which encompasses, in particular, 
the three integer values d — 2, 3, and 4 (compare Fig.|3|). Thus the curvature parameters 
are overdetermined. In the present case, however, this consistency may be regarded as a 
simple consequence of the differential equation ( |4.3| ) in which the perturbation generated 
by the surface curvature is proportional to d — 1. In the next subsection we consider the 
small curvature expansion (|3.12j ) for the excess adsorption, for which the corresponding 
consistency check in D = 4 provides a more stringent test. 

(iii) y± — > with s and £± fixed, i.e., R — > (compare Sec. pCp . We consider the case 
of a sphere, i.e., d = D with D = 4. By inserting the ansatz m(s; R, r) — m± i j ) = Ru(s; r) 
into Eq. (|4.3p and keeping only terms linear in R one obtains a linear differential equation 
for u(s; r), similar as in case (i) above. The solution of this differential equation in terms 



of the modified Bessel function K\{x±) in conjunction with Eq. (17?) leads to 



P ± (x ±} y ± ^0) - P ±>b -> 2c±y ± x± 1 K 1 (x ± ) , d = D = 4 . (4.18) 



Here the constant prefactor 2c± is fixed because the limit x± — > of the rhs of Eq. (|4.18 



must reproduce Eq. ( |2.12j ), in which C±(j — > 00) — > 2c±/^y for d = D = 4 according to 
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Eq. ( |2.15D . Equation ( |4.18| ) is in line with Eq. ( |2.31| ) because for D = 4 one has (3/v = 1 



and the universal scaling functions J 7 ±(x±) are given by 

F±(x±) = x±K x {x±) , D = A. (4.19) 

B. Excess adsorption 

For D = 4 Eq. ( |3.7| ) reduces to 

T(t^0,R) -> Aa^{g ± \ln\t\ | + G±(y±)} , D = 4 , (4.20) 



with g + = l/v2 and g_ = 1. The universal scaling function G±(?/±) in Eq. (|4.20 ) is given 



by Eq. ( |3.6|) with P±(x±,y±) and P±(x±) from Eqs. (|4.7|) and (|2.9|) . In the following we 
discuss the behavior of G±(y±) in the limit of large and small values of y± = R/£± and 
give results for the whole range of y± representing the crossover between these two limits. 

(i) G±(y±) for y± — > oo (compare Sec. |III B| ) . By using the differential equation (|4.3|) 
in conjunction with Eq. ( |4.7| ) one can determine the coefficients a^ D and bf D in Eq. (|3.11|) 
for generalized cylinders K with D = 4 and arbitrary d in the interval 1 < d < 4. This 
calculation is presented in Appendix |C[ As a result one finds that the dependences of a~f A 
and on d are precisely of the form given by Eq. (|3.13|) . In the present case D = 4 one 
has Itl 13 ' 1 ' = 1 and the curvature parameters are given by [[16[ 



A+ = 5.09 , Af = 4.91 £_ , (4.21a) 



A+ = -1.55 , A 2 - = -0.91 e , (4.21b) 

A^ = 2.87 , \ G = 2.52 ^ 2 . (4.21c) 

The consistency with Eq. ( |3.13| ) for D = 4 can be traced back to nontrivial properties 
of the differential equation ( |4.3| ) (see Appendix 0) and thus provides a valuable and 
important check for the validity of the general curvature expansion ( |3.12| ). We expect 
that the curvature expansion ( p. 12] ) is also valid in D = 3 and that the corresponding 
numerical prefactors multiplying £ or £ 2 in the curvature parameters for D = 3 differ only 
quantitatively from those in Eq. ( |4.21|) valid for D = 4. 
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As an illustration, consider a curved membrane with both sides exposed to a fluid near 
criticality. In particular we consider the total excess adsorption, i.e., the sum of the excess 
adsorptions on each side of the membrane, per unit area. In this case the contributions 
to X^K m in the expansion (|3.12|) from each side cancel and the signs in Eqs. ( |4.21b|) and 
( |4.21c ) in conjunction with Eq. ( 2.21 ) imply that the total excess adsorption is larger 



near spherical regions (d = D) of the membrane as compared to flat regions, whereas 
near cylindrical regions (d = 2) it is smaller as compared to flat regions. 

(ii) G±(y±) for y± — > (compare Sec. p^IIC]) . In this case the behaviors for spheres 



and cylinders are qualitatively different due to the different behaviors of P±(x±, y±). For 
a sphere (i.e., d = D with D — 4) the power law ( |3.16| ) is valid where the exponent 
— D + 1 + P/u equals —2 and according to Eqs. ( |3.17| ) and fl4.19| ) the universal amplitudes 
uj± are given by 



A\/2 , u_ = 8 , d = D = A . (4.22) 



For a cylinder (i.e., d = 2 and D = 4) one finds the behavior ( |3.20|) where the universal 



numbers v± in Eq. ( |3.21 ) can be evaluated numerically with the result 



v + = 1.90 , t>_ = 1.86 , d = 2, D = 4 . (4.23) 

(iii) The full scaling functions G±(y±) describe the crossover between their analytic 
behaviors for y± = R/£± — > oo and the power laws for y± — > which have been discussed 
in (i) and (ii) above. Figures |8| and |S] show numerical results for G±(y±) corresponding to 
a sphere and a cylinder, respectively, in D = 4. The results confirm the small curvature 
expansion as implied by Eqs. (|3.11|) - Q3.13Q , and (^4.21|) and the properties ( |3.16|) and 



( |3.20| ), and they provide the range of validity of the asymptotic behavior. Note that in 
the case of the sphere y±G±(y±) as function of y^ 1 = £±/R diverges for y^ 1 — > oo with 
the corresponding power law whereas in the case of the cylinder with increasing y^ 1 it 
interpolates between one finite value related to X x to the other finite value v±. 



FIG. 8. Scaling functions G±(y±) for a sphere in mean-field approximation (i.e., d = D with 
D = 4) as function of y^ 1 = £±/R for (a) T > T c and (b) T < T c . The dotted lines show 
the small curvature expansion (see Eqs. ( pl|) - ( gig ), and ( gjj) )) valid for y ± 1 < 1 and the 
dashed lines show the power law (see Eqs. ( |3.16 ) and ( 4.22[) ) valid for y^ 1 — > oo; for D = 4 the 
exponent —D + 2 + j3/v is equal to —1. 



34 




35 




FIG. 9. Scaling functions G±{y±) for a cylinder in mean-field approximation (i.e., d = 2 and 
D = 4) as function of y^ 1 = £±/R for (a) T > T c and (b) T < T c (compare Fig.||). The dotted 
lines show the small curvature expansion and the dashed lines correspond to the numbers v± in 
the behavior G±{y± -> 0) -> v± yg 1 (see Eqs. ( |3T2C| ) and flE23p). 
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V. SUMMARY AND CONCLUDING REMARKS 



We have studied critical adsorption phenomena on spherical and cylindrical particles 
of radius R which are immersed in a fluid near criticality, t = (T — T c )/T c — ► 0, for the 
case that the fluid is at the critical composition. The corresponding adsorption profiles 
at the radial distance s from the surface are characterized by universal scaling functions 
P±( s /£±, R/£±) f° r T 7^ T c , involving the bulk correlation lengths £± for T ^ T c , and 
C + (s/R) for T = T c (see Eqs. ( f2.1| ) and ( P-llQ , respectively). 



In the following we summarize our main results starting with local properties of 
P±{x±,y±) an d in various limiting cases as indicated by Fig.f|. 

(1) For T = T c we have introduced the short distance expansion of the order parameter 
profile near a weakly curved (D — l)-dimensional surface of general shape (see 
Eq. ( |2.18| )). This expansion involves the local curvature invariants K m , K^, and Kq 
(see Eq. (|2.17Q ). The corresponding expansion parameters Ki, K2, and appear 
also in the short distance expansion (|2.23| ) of P±(x±,y±) valid for T ^ T c . The 
parameters k±, k 2 , and kq are universal and depend only on the space dimension D 
(see Eqs. flD3D , ( ggg ), and fliTTD ). 



(2) For R -C s, £ the order parameter profile near a cylinder becomes independent 
of i? in the limit R — > 0, i.e., if the cylinder radius is microscopically small (see 
Eqs. ( p.33| ) - ( |2.36| )). In contrast, near a sphere the 'universal part' of the order 
parameter profile, i.e., the one described by Eqs. ( |2.30p and ( |2.31| ), vanishes for 
R — > 0. This qualitative difference in behavior can be explained by means of a 
small radius operator expansion (see Eqs. (|2.24j ) and ( |2.25|) ). As indicated by Fig. |] 
a sphere is an irrelevant perturbation for the fluid near criticality whereas a cylinder 
is a relevant perturbation. 

(3) The explicit forms of the universal scaling functions P±(x±, y±) and £+(7) for a 
sphere and a cylinder within mean- field approximation (see Figs. |5| - 0) corroborate 
the results described in the above paragraphs (1) and (2). 



We now turn to the excess adsorption T(t —>■ 0,R) describing the total enrichment of 
the preferred component of the fluid near criticality in proximity of a sphere or cylinder 
(see Eq. (|3.1|) ). The curvature dependence of T{t — > 0,R) is characterized by universal 
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scaling functions G±(R/£±) obtained from P±(s/£±, R/£±) by integrating over the first 
variable s/£± (see Eqs. ( |3.6| ) and (|3.7|)). 



(4) For R/£± > 1 we have introduced the expansion of T(t — > 0, R) near a weakly 
curved (D — l)-dimensional surface of general shape in terms of the local curvature 
invariants K m , K^, and K G (see Eq. ( p,12|) ). The corresponding expansion param- 
eters Af 1 , A^, and Xq depend only on the space dimension D and can be expressed 
in terms of the universal coefficients a^ D and b~f D appearing in the expansion of 
G±(y± — > oo) (see Eqs. ( |3.11| ) - ( |3.13j )). The explicit calculation of A-jf , A^, and 



A^ within mean-field approximation (see Eq. (4.21) and Appendix |C]) provides an 



important check for the validity of the general curvature expansion of the excess 
adsorption. 

(5) For R/£± — > the behaviors of T(t — > 0, R) and G±(y±) are characterized by power 
laws (see Eqs. (|3.16|) - ( 3.18|) for a sphere and ( p.20|) - Q3.22j ) for a cylinder). The 



exponents in these power laws are different for spheres and cylinders beyond the 
purely geometrical effects. This can be traced back to the corresponding difference 
in behavior of P±(x±, y±) addressed in paragraph (2). 

(6) The explicit forms of G±{y±) for a sphere and a cylinder in mean-field approximation 
(see Figs.|| and ^) confirm the general results presented in paragraphs (4) and (5). 

We conclude by summarizing some of the fieldtheoretic developments and perspectives 
for future theoretical work, simulations, and experiments. 

(7) For a systematic field-theoretical analysis of the critical adsorption phenomena on 
spheres and cylinders it is useful to introduce the particle shape of a 'generalized 
cylinder' (see Eq. ( |1 . 1| ) and Fig. |1]) which is characterized by the space dimension D 
and an internal dimension d encompassing a sphere, a cylinder, and a planar wall 
as special cases. 

(8) For a single sphere, certain asymptotic behaviors of the order parameter profile are 
characterized by universal quantities (see Eqs. ( |2.23| ) and ( |2.31| )). These quantities 



are accurately known in D = 3 (see Eq. ( |2.19p , (2.20), and Table ^ as far as the 



scaling functions J r ±(x±) are concerned, only the function jF + (s + ) appropriate for 
T > T c is known accurately for D — 3, see Appendix 0). The corresponding 
predictions can be tested using small angle scattering of light, X-rays or neutrons 
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which can probe the enlarged effective size of the colloidal particles due to the 
adsorption. 



(9) For a single sphere, in D = 3 the behavior of the excess adsorption in the limit 
R/i± — ^ is available as well (see Eqs. ( |3.7|) , ( |3.16|) , and Table p|) and can be 



tested experimentally by, e.g., volumetric measurements. For a disc in D = 2 the 
numerical values of the universal amplitudes u± in Eq. ( |3.16| ) are also known (com- 
pare Appendix Bf). This can be relevant for protein inclusions in fluid membranes 



(10) Apart from the curvature dependence of the excess adsorption (see Eq. ([3.12 )) also 



the corresponding energetic contributions, i.e., the change of the surface tension 
generated by the surface curvature, can be relevant for applications. For example, in 
case of a membrane immersed in a fluid near criticality such contributions and their 
temperature dependence are expected to influence the intrinsic bending rigidities of 
the membrane. These modifications of the bending rigidities can, in turn, induce 
shape changes of vesicles formed by closed membranes in a controllable way (see, 
e.g., Ref. HI). 



11) A one-dimensional extended perturbation in an Ising-like system which breaks the 
symmetry of the order parameter represents a relevant perturbation (see Sec. [II C ) 



of the bulk system. This gives rise to new universal quantities such as the criti- 
cal exponents m and rj^ characterizing the decay of the structure factor (see Ap- 
pendix [A] and in particular Eqs. (|A9Q and ( |A11| )) and the amplitudes n± and v± (see 



Eqs. ( [2.360 an d ( |3.21| ), respectively; the mean- field values of n± and v± are quoted 
in Eqs. ( (4.10[ ) and ( |4.23| )). Apart from a rodlike particle immersed in a fluid, such 



a one- dimensional perturbation could also be realized in a solid, e.g., by a lattice 
dislocation in a binary alloy if the dislocation locally breaks the symmetry of the 
order parameter. This would be a natural starting point for testing the numerous 
fieldtheoretic predictions derived here by Monte Carlo simulation. 

^12) As stated in the Introduction, the results addressed in paragraph (2) can be rel- 
evant for the flocculation of colloidal particles dissolved in a fluid near criticality. 
Specificly, if one considers the case T = T c and compares the behaviors of the order 
parameter profiles near a sphere and near an infinitely elongated cylinder, respec- 
tively, the order parameter of the fluid at a distance s ^> R from the particle with 
radius R is larger by a factor ~ {s/R) 13 ^ 1 for the cylinder than for the sphere 
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(see Eqs. ( |2.11| ), ( |2.30| ) and ( |2.33| )). For rods with a large but finite length I and 
distances s > I this enhancement is of the order of (l/R) l3 l u ^> 1. Since the critical 
Casimir forces between particles are partially induced by an overlap of the order 
parameter profiles generated by the individual particles in the space between them, 
the aforementioned results suggest that not only spheres but, a fortiori, also long 
rods in a fluid near criticality can aggregate due to critical Casimir forces. 
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APPENDIX A: TWO-POINT CORRELATION FUNCTION NEAR A NEEDLE 

The two-point correlation function G(r, r ') for Gaussian fluctuations around the profile 
m(r±) = l/r_|_ (see the text following Eq. ( [4. 1 3| ) ) at criticality satisfies the differential 
equation 

- A D + 3 m 2 (r ± )j G(r, r ') = 5 (D) (r - r ') (Al) 

where is the Laplacian operator in Z)-dimensional space. (Although the present mean- 
field calculation implies D = 4 we leave the symbol D for clarity.) In coordinates adapted 
to the geometry of a generalized cylinder K with d = 2 and R = the Laplacian operator 
Ajj has the form 

D-2 ~ 2 

A D = A ± + qTY (A2) 

i=l II. i 



with 



d 2 1 d 1 d 2 

+ - IT" + Z2 IZZ (A3) 



dr]_ rj_ <9r_|_ r]_ dd 2 

as the Laplacian operator in the two-dimensional subspace perpendicular to K. Here d 
denotes the angle between and a fixed direction in the radial subspace. In order to 
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I 




FIG. 10. A thin 'needle' corresponding to a generalized cylinder K with d = 2 and R = 
(compare Eq. and Fig. ([!])). The two spatial arguments r and r' of the correlation function 
in Eq. (KA) are also shown. 



solve Eq. (|A1|) we carry out the Fourier transform in the subspace parallel to K and a 
partial wave decomposition in the radial subspace, i.e., 

G(r,r') = G(r ± , r;,tf, I'll I) (A4) 
= W ™(®) I (2tt\d-2 ex P[^ P ( r ll - r || )] G n {r ± ,r±,p) , 

n=0 ^ 71 ' 

where i? is now the angle between r_|_ and (see Fig. 110). The functions W n {$) are given 
by 

W n {f) = 2 ~ * ra '° cos(ntf) (A5) 
with <5 n)0 = 1 for n = and zero otherwise. They obey 'Yl n W n ('d) = 8{&)/2 where 
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the support of the 5-function on the rhs is located entirely in the interval $ > so that 
£ d$5($) = lfor any d 1 > 0. The propagator G n in Eq. (|A4]) satisfies the radial equation 



d 2 Id 2 n 2 + 3 
+ V + 



dr 2 rj_ dr± 1 r 



,2 



Gn(rx,r±,p) = . (A6) 



For r± ^ r[_ the solutions of this equation are modified Bessel functions K\(pr±) and 
1 \{p r±) |^5] with index 

A = A n = Vn 2 + 3 . (A7) 



The physically acceptable solution G n of Eq. ( |A6[ ) is uniquely determined by the require- 
ments that (i) it is symmetric in r± and r ± , (ii) it decays for r± — ► oo, and (iii) it does 
not diverge as the singularity at r± = in Eq. (|A6|) is approached. The result is 



G n (r ± ,r ± ,p) = K x (prP)I x (prf) (AS 



with = max(r_|_,r^ ), r± = min(r_|_,r^ ), and A given by Eq. (|A7|) . Equations (|A4|) 
and ( |A8|) represent the two-point correlation function near a thin needle in the presence of 
the adsorption profile m(r±) = l/r± at criticality. In the following we infer the asymptotic 
behavior of G(r,r') in various limits. 



(a) r±, r ± fixed and l r || — r || I — > oo- Therhsof Eq. ( |A5| ) can be expanded for p — > 0, e.g., 
by relating K\ to J_a and /a via the corresponding formula in section 9.6.2 of Ref. [45(a)] 
and by using the formula in section 9.6.10 of the same reference. Apart from terms pro- 
portional to integer powers of p 2 , which are analytic in p yielding only exponentially 
decaying contributions on the lhs of Eq. ( |A4| ) , one finds that the leading singular contri- 
bution behaves as p 2X . Regarding the leading behavior for |r|| — r,, | — ► oo in Eq. (|A4|) only 
the term with n = is important so that 



G(r ± , r{_ , tf, | r n - \ oo) ~ _ - -. D _ 2 + m ( A9a ) 



with 

^ = 2A = 2\/3 . (A9b) 



(b) r || = r,i , r ± fixed and r± — > oo. For r± > r ± using Eq. (|A4|) and rescaling x = pr± 
leads to 
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1 00 f 



d " K x (x)I x (xrJr ± ) . (A10) 



(2tt 



D-2 



Here I\(xr±/r±) can be expanded for r±/r± — > so that in this limit again only the term 
with n = is important and 

G(r ± -> oo, r; , tf, 0) ~ ^4^7 (Alia) 

r ± 

with 

77-l = A = V3 . (Allb) 

We note that within the present mean-field theory the relation 2r]± — T] + m is fulfilled 
because T](D = 4) = 0. 

The exponents in Eqs. (|A9|) and ( |A11|) for a thin needle are different from those for 



the half-space for which r]^ s = 6 and ^ s = 3 in mean-field approximation. Thus the 
correlations in the critical fluid are more suppressed in proximity of the needle than in 
the bulk fluid, and are even more suppressed near the surface bounding the half-space. 
This can be understood by noting that a needle represents a weaker, but still relevant, 
perturbation for the critical fluid than a planar wall (compare Fig. ^). For a sphere one 
finds that the two-point correlation function at criticality decays as |r — r / |~( D ~ 2 +'?) with 
the bulk exponent r] if the distance vector r — r ' is increased in any direction. This reflects 
the fact that a sphere represents an irrelevant perturbation. 



APPENDIX B: AMPLITUDES u± AND SCALING FUNCTIONS F±{x±) 

In this appendix we determine the universal amplitudes u± (see Eqs. ( [3.1 6|) and ( |3.17| )) 
and the universal scaling functions J-±(x±) (see Eqs. (|2.31|) and Q2.32Q ). The corresponding 
mean-field results (see Eqs. (|4.22|) and ( |4.19|) , respectively) hold in D = 4. For D = 3 and 
2 results are available p3|-p)4| for the Fourier transform of the bulk two-point correlation 
function introduced in Eq. (|2.32[) , i.e., 



Stip) = J cPre'*'<*(r)<&(0))S , (Bl) 

which can be written in the scaling form 

S t (p) = ^{trg^p^) . (B2) 
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The universal scaling functions g±(Y±) are fixed by the normalization condition g±(0) = 1 
and by the choice of £± as the true correlation length for T ^ T c . The amplitudes C are 
nonuniversal and 7 = v{2 — rj) is a universal bulk critical exponent. Numerical values of 
the universal bulk critical exponent rj are given by rj(D = 2) = 1/4 and rj(D = 3) = 0.031 
p3| (compare Eq. (|2.5|) ). Inserting Eq. ( |2.32| ) into Eq. ( [ETD and using Eq. ( |3.17| ) yields 

u± = (B3) 
QdQ± 

with flu = 2ti d / 2 / T(D/2). The universal numbers Q± are given by 

N Rl~ v Q 3 , T>T C , 

Q± < r+ (B4) 

NR 2 + ^Rf 2 ^Q 3 , T<T C , 

with the universal amplitude ratio Q3 = D (£o + i) 2 ~VC + introduced by Tarko and Fisher 
[ 51| . The amplitudes £q\ correspond to correlation length defined via the second moment 



HJ of the correlation function, i.e., £±,i = ^cnl^l"^- The amplitude D is defined by 
St=o(p) = Dp 71 ' 2 which implies that Bq> in Eq. ( |2.27|) is related to D via Bq> = ND with 
the numbers 

r(|)/[2vr2 3 / 4 r(|)] = 0.654308 , D = 2 , 
N = { N (B5) 

r(?7) sin(^)/(27r 2 ) = 0.078196 , D = 3 . 

Numerical values of Q 3 are given by Q 3 (D = 2) = 0.41377 and Q 3 (D = 3) = 0.896 
51| . In Eq. ( |B4|) R + = £q /C^, R$ = ^o"/^o", and C + /C~ are known universal amplitude 



ratios with R + (D = 2) = 1.000402 (see Eq. (3.7) in Ref. 0), R + (D = 3) = 1.0003 |5( 
Rt(2) = 2, i? ? (3) ~ 1.92 0; (C r+ /C r -)(2) = 37.693562 0, (C + /C _ )(3) = 4.95 ^ 
Using Eqs. ( |B3[ ) and ( P^ ) in conjunction with the values of c±(D = 2) from Table [I] and 
with c + (D = 3) = 0.94, and C-(D = 3) = 1.24 yields the values of uo± quoted in Table 



III. We note that the accuracy of the quoted value of uj-(D = 3) is unknown because the 



accuracy of the aforementioned value of R^(D = 3) is not given reliably pD 



Next we outline how the universal scaling functions J r ±{x±) introduced in Eq. ( |2.32| ) 
can be inferred from presently available results for g±(Y±) (see Eqs. ( pl| ) and (P2j) ). 

(a) D = 3: In this case the scaling functions F±{x±) are given by 

00 

F±(x±) =k±xl J dY± Y± sm(Y ± x ± ) g ± (Y±) , D = 3 , (B6) 



44 



with the amplitudes k± fixed by the condition .F±(0) = 1 which allows one to express the 
nonuniversal amplitudes in terms of B§ and £q • F° r the case T > T c the approximation 
for g + (Y + ) proposed by Bray j53f can be regarded to be reliable. Accordingly, for values 
Y+ — 20 the function g + (Y + ) can be inferred from the first column in Table V of Ref. [|53|1 
whereas for Y + > 20 the asymptotic expansion by Fisher and Langer [J5| is applicable, 
i.e., 



g + (Y + -> oo) 



2—77 



1 + 



Co 



Y 



(l-a)/i/ 



+ 



v 



1/1/ 



+ 



7J 



(B7) 



with coefficients Cj and the bulk critical exponents rj, a, and v. Using in Eqs. (|B6|) and 
( p7|) the values 77 = 0.041, 1/ = 0.638, and a = 0.086 as quoted in Ref. |53| yields 
C l = 0.909, C 2 = 3.593, C 3 = -4.493 p[, and £;+ = 0.7166. Unfortunately, for T < T c 



and .D = 3 we are not aware of an accurate estimate of g-{YJ). 

(b) D = 2: In this case the scaling functions J r ±(x±) are given by 

00 

T ± {x ± ) =k' ± xl J dY ± Y ± J {Y ± x ± )g±{Y ± ) , D = 2 



where Jo is a modified Bessel function [[45 
deduced from Ref. 



In D = 2 exact results for g±(Y±) can be 
. However, some care is necessary regarding the definition of the 
scaling variable y used in Ref. [p2[j: whereas y = Y + for T > T c , one has y = 1.959 Y- 

1 with S2 from 



for T < T c . The numerical prefactor in front of F_ equals (yE 2 i?_ 



Table IV in Ref. |52| (see Eq. (3.7) in Ref. ||52|| ) and with the universal amplitude ratio 
i?_ = £o7£(Ti = 1-615 in D = 2 ]5T| . For the amplitudes fc ± in Eq. ( P8| ) one obtains 
k' = 0.5874 and k' = 0.05055. 



APPENDIX C: SMALL CURVATURE EXPANSION OF THE EXCESS 

ADSORPTION 



In this appendix we outline how the coefficients a dD and b dD introduced in Eq. (|3.11|) 



can be calculated within mean- field theory, i.e., for D = 4. Since the approaches for 
T > T c (+) and for T < T c (— ) are quite similar we restrict our presentation to the case 
T > T,. 



According to Eqs. (|4.3|) and ( f4.7|) for D = 4 the universal scaling function P + (x + ,y + ) 
satisfies the nonlinear differential equation (in the following we drop the subscript '+' in 
order to simplify the notation) 
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d-1 



P"{x,y) + — P'(x,y) - P(x,y) = P 3 (x,y) 



(Cla) 



x + y 

where the derivatives are taken with respect to the variable x. The parameter d with 



1 < d < 4 can be chosen arbitrarily. Equation (|Cla ) is supplemented by the boundary 
conditions 

P(x^0,y) -> ^- , P(oo,y) = P + , b = . 
x 

The scaling function = G+(y+) from Eq. ( |3.6| ) then reads 



(Clb) 



G(y) = j rfx{(- + l)^ 1 p(x,y) - P (x)} 
o 



(C2) 



with P(x,y) defined by Eq. (|C1]) and the half-space profile Po(o;) = P(x, oo) = P + (a; + ) 
given by the first part of Eq. ( |2.9| ). Note that the singularity at x = of the first term 
in curly brackets in Eq. ( |U2| ) is cancelled by the second term. According to Eq. (3.11) the 
function G(y) can be expanded as 



G(y -> oo) = a d y 1 + b d y 2 + ... 



(C3) 



with the coefficients = and 6^ = 6/4 which we want to determine. To this end we 
expand P(x, y) as 



P(x,y^ 00) = P (x) + P^x)?/" 1 + P 2 (x)y + 



(C4) 



By inserting Eq. ( |U4| ) into Eq. ( |Cla| ) and equating terms with the same power in y one 
derives the familiar nonlinear differential equation for the half-space profile Pq(x), i.e., 



tjH r> p3 

mi ~ -ni - -mi 1 



Po(x -> 0) 



v/2 



Pn(oo) = 



(C5a) 



(C5b) 



and the two linear differential equations 

P» - Pl - w 2 P 1 = -(d-l)P^ , 

Pi(0) = ~{d-l) , Ptioo) = 
o 



(C6a) 
(C6b) 



46 



and 



P'i - P 2 - 3P 2 P 2 = -{d - l)P[ + 3P P 2 + (d-l) xP^ , (C7a) 

^(0) = - || (d - I) 2 + ^ (d - 1) , ^2(00) = , (C7b) 

for P\{x) and P 2 {x), respectively. For the boundary conditions in the second parts of the 
above equations compare Eqs. ( [4.1 6| ) and ([4.7|). By using 



(- + l) d 1 = 1 + (d-l)xy- 1 + l(d-l)(d-2)x 2 y- 2 + ... (Q 



in the integrand on the rhs of Eq. ( p2|) in conjunction with Eq. ( |U4j) we find for the 
coefficients in Eq. ( |C3| ) the expressions 

00 00 



and 

00 

6 d = \){d-2) J dxx 2 P (x) (CIO) 



00 00 


where PqC^), Pi(^) 3 an d -^(^c) are the solutions of Eqs. QC5| ) - (|C7[). While Po(^) is 
given by the first part of Eq. (|2.9| ) we now turn to the calculation of P\(x) and P2(x). 
It is important to retain the dependence of the latter functions on d in analytical form 
in order to be able to carry out the consistency check of Eq. ( |3.13| ) for the present case 
D = 4. 

(i) Function P\(x). Both the rhs of Eq. ( p6aj) and the boundary condition in the first 
part of Eq. ( |C6b|) depend on d via the term (d — 1). We note that even the full function 
Pi(x) exhibits this simple dependence on d. The reason is that Eq. ( |C6a| ) is linear with 
respect to Pi so that Eq. (03) is solved by the ansatz 



p x (x) = (d-l)U(x) (Cll) 
where U(x) is independent of d and satisfies the inhomogeneous linear differential equation 
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(C12a) 



*<0) = 



U(oo) = . 



(C12b) 



The solution of Eq. ( |C12| ) can be given analytically (see, e.g., Ref. 
cosh(x) — sinh(a;) 



U(x) 



x 



6v^2 sinh 2 (a;) 

2 cosh(2x) — 2 — 3x — 3x cosh(2x) + 3 sinh(2x) — 3x sinh(2x) 



(C13) 



(ii) Function P%{x). By using Eq. ( |C11| ) the differential equation ( |C7a| ) turns into 
P'i - P 2 - 3P Q 2 P 2 = -(d - I) 2 U' + (d - if 3P U 2 + (d - 1) xPq (C14) 



with U(x) from Eq. (|C13| ). In this case both on the rhs of Eq. (|C14| ) and in the boundary 
condition in the first part of Eq. (|C7b|) different dependences on d arise via the terms 
(d — l) 2 and (d — 1). Again it is crucial to observe that Eq. (|C14 ) is linear with respect 
to P 2 so that its dependence on d takes the simple form 



P 2 (x) = (d-l) 2 V{x) + (d-l)W(x) 



(C15) 



where V(x) and W(x) are independent of d and satisfy separately the inhomogeneous 
linear differential equations 



V" - V - 3P 2 V 



-U' + 3P n U 2 



(C16a) 



and 



^(0) = -|f, V{oo) 



W" -W - 3P 2 W = xPq 



W'(0) 



(C16b) 



^ , W(oo) 



. 



(C17a) 



(C17b) 



Equations (|C16|) and (|C17|) can be solved numerically, which is facilitated by the fact 
that both U(x) and Po(x) on the rhs of Eqs. (|C16a| ) and (|C17a ) are known in analytical 
form. By inserting the resulting function P^ix) from Eq. ( |U15 ) and the function P\{x) 
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from Eq. ( |C11| ) into Eqs. ( |U9| ) and ( |C10| ) one obtains the dependences of a^~ 4 = and 
b£ 4 = ha on d as given by Eq. ( |3.13| ) with the curvature parameters given in the first parts 
of Eq. (|42lD . Theproced ure in the case T < T c is completely analogous to that for T > T c 
as outlined above and yields the curvature parameters in the second parts of Eq. (|4.21| ). 
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